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We develop a model for the effect of a magnetic field on quasiparticle interference in an iron- 
based superconductor. Recently, scanning tunneling experiments have been performed on Fe(Se,Te) 
to determine the relative sign of the superconducting gap from the magnetic-field dependence of 
quasiparticle scattering amplitudes. Using a simple two-band BCS model, we study three different 
cases of scattering in a spin-split spectrum. The dominant effect of a magnetic field in iron-based 
superconductors is caused by the Pauli limiting of conduction electrons. Thereby time reversal odd 
scattering is induced which enhances the sign-preserving and depresses the sign-reversing peaks in 
the quasiparticle interference patterns. 
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I. INTRODUCTION 

The discovery of iron-based layered pnictide and 
chalcogenide superconductors^ with superconducting 
(SC) transition temperature as high as 55K (Ref.— ) has 
generated enormous interest in the physics of these ma- 
terials. Similar to the cuprates, the iron-based supercon- 
ductors are highly two-dimensional and superconductiv- 
ity occurs in close proximity to anti-ferromagnetic order^ 
leading several groups to propose that the pairing is 
driven by anti-ferromagnetic spin fluctuations 

The first step in identifying the pairing mechanism is 
to investigate the structure of the SC-gap function, which 
describes the strength and quantum mechanical phase of 
electron pairs in momentum (k) space. While the SC 
gap function of conventional phonon-mcdiatcd supercon- 
ductors has the same phase throughout momentum space 
(s-wave symmetry) , that of spin- fluctuation mediated su- 
perconductors is expected to exhibit a sign reversal be- 
tween those Fermi momenta connected by characteris- 
tic wave vector Q of the spin fluctuations £ In d wave 
superconductors, the nodal planes of the order parame- 
ter intersect with the Fermi surface, leading to gaplcss 
quasiparticle (QP) excitations that can be detected ther- 
modynamically and by low energy probes. However, if 
the sign reversal occurs between disconnected hole and 
electron pockets (compare Fig. Q]) resulting in an "s± 
symmetry,"— £ the relative sign of SC gap must be de- 
termined by phase-sensitive experiments such as Joseph- 
son junctions experiments^ or composite SC loopsJ^ An 
alternative technique is scanning tunneling microscopy 
(STM) which determines the dispersion of quasiparticle 
states from the quasiparticle interference (QPI) patterns 
induced by impurity scatterin g 11 15 . Such experiments, 
performed in an external magnetic field, offer the capa- 
bility of probing the phase of the superconducting order 
parameter by detecting a field enhancement of the sign- 
preserving scattering that results from the sensitivity of 
QPI to the coherence factors associated with impurity 
scattering. 



The recent detection of a field enhancement of quasi- 
particle scattering between the hole pockets of the 
iron chalcogenide superconductor Fe(Se,Te) has been 
interpreted^ as evidence for the s± pairing scenario. 
However, there are certain differences between the field 
effect in iron-based and cuprate superconductors that 
raise questions about this interpretation. For exam- 
ple, Fe(Se,Te) is a Pauli-limited superconductor where 
the coupling of the field to the conduction electrons be- 
comes important. In this material, the Zeeman split- 
ting is a large fraction of the superconducting gap size, 
firiB/A « 0.4, generating new components to the field- 
induced scattering that have been hitherto neglected in 
models of quasiparticle interference. 

In the cuprates, quasiparticle scattering in external 
magnetic fields is strongly affected by vortices. How- 
ever, STM measurements on Fe(Se,Te) did not observe 
appreciable correlation in space between the location of 
vortices and the magnitude of field-induced change in 
the QPI intensity^ These results suggest that in Pauli- 
limited superconductors the role of vortices is much less 
important compared to homogeneously distributed field- 
induced scatterers created by Zeeman splitting. One of 
the key arguments advanced in the STM measurements 
of Fe(Se,Te), is that an applied magnetic field will, in 
general, enhance the time-reverse odd components of the 
scattering that lead to sign-preserving scattering. How 
effective are these arguments in the presence of signifi- 
cant Pauli limiting effects, and to what extent can we 
attribute the observed field enhancement of the sign- 
preserving scattering to the Zeeman splitting? 

Motivated by the above considerations, this paper de- 
velops a phenomenological model for the ficld-dcpcndcnt 
QPI in iron-based superconductors, taking into account 
the Zeeman splitting of the quasiparticle dispersion. Us- 
ing a simple two-band BCS model for an unconventional 
s-wave superconductor, we calculate the STM conduc- 
tance ratio for different types of disorder in an external 
magnetic field at which we neglect the scattering off vor- 
tices. In our model, we assume an s±-wave symmetry for 
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the superconducting order parameter which means that 
the size of the superconducting gap is isotropic along the 
hole and electron pockets. Our study confirms that Zee- 
man splitting provides a strong source of time-reversal- 
odd, sign-preserving scattering. We also find that if res- 
onant or magnetic scatterers are present in the super- 
conducting material, the spin split spectrum leads to an 
enhancement of sign preserving scattering. QPI patterns 
for non-magnetic and magnetic impurities without Zee- 
man field have been computed theoretically by a number 
of works using more realistic band structures! 17 21 

The paper is organized as follows. In Sec. II A, we 
briefly resume the basic concepts behind some types of 
scatterers with different time reversal symmetry proper- 
ties and their coherence factors. The BCS model for an 
iron-based superconductor will be described in Sec. II B. 
After a short discussion about calculating the tunneling 
conductance from local density of states in Sec. Ill, we 
present our results of the numerical evaluation in Sec. IV. 
We conclude in Sec. V. 



We now review the connection between coherence fac- 
tors for scattering in a BCS superconductor and the time- 
reversal symmetry of the scattering process. Consider a 
scattering process described by the i-matrix, 
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where the Nambu-spinors Vt^ and 'Iv account for differ- 
ent particle states according to Eq. (fTJ) . In the supercon- 
ducting state, the i-matrix determines the Friedel oscil- 
lations in the local density of states that are probed in an 
STM experiment. To determine the matrix elements for 
quasiparticle scattering, we apply a Bogoliubov transfor- 
mation to the particle spinors in Eq. ([3]), leading to 
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II. PHENOMENOLOGICAL MODEL 
A. Coherence factors 

An STM measurement probes quasiparticle interfer- 
ence (QPI) by observing Friedel oscillations in the tun- 
neling density of states induced by impurity scattering. 
The Fourier transform of these fluctuations in the density 
of states provides information about the rate at which 
quasiparticles are scattered by impurities. In general, 
the rate of impurity scattering is proportional to the co- 
herence factors associated with the underlying scattering 
mechanism. By measuring the momentum dependence of 
the scattering, it becomes possible to extract phase sen- 
sitive information about the underlying order parameter. 

Quasiparticles in a superconductor are a coherent su- 
perposition of electrons and holes. Coherence factors de- 
scribe the difference between the scattering rates of a SC 
quasiparticle and a bare electron off a given scattercr^ 
In general, the scattering rate of an electron from ini- 
tial momentum k to final momentum k' is determined 
by an energy-dependent t- matrix T(k, k',a>). The t- 
matrix, which depends on the quasiparticle spin degrees 
of freedom is most conveniently described using four- 
component Nambu notatio n 23 ; 24 . For a single impurity 
with scattering potential ?7(k, k') in momentum space, 
the corresponding i-matrix exactly accounts for multiple 
scattering off that impurit y 24 ) 25 In first-order perturba- 
tion theory (Born approximation) the i-matrix is equal 
to U (k, k'). In four-component Nambu notation the elec- 
tron field is described by a four-component spinor con- 
taining particle and hole operators^ 3 - 
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where for convenience, we have assumed a gauge in which 
the Bogoliubov coefficients are real. For a BCS super- 
conductor with electronic dispersion £k, gap function Ak 
and quasiparticle energies 
ubov coefficients arc 



+ A 2 , the Bogoli- 




^ = sign(A k )Ji h-g 



(4) 
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Introducing a Nambu spinor for the quasiparticle op- 
erators <i> k = (a k a^ k ^ a- k .| — a k,t), we may write the 
Bogoliubov transformation ([3]) as 



*k = U«k 1 + i Vk fa) $ k , 



(0) 



where 1 is the 4x4 unit matrix and To, is the isospin 
matrix, 



To = / 



-1 

1 



(7) 



where an underscore denotes a two dimensional matrix. 
Using Eqs. © and ^ the i-matrix can be written in the 



form T = E k k , f££ $k' , where 



fQ p - 

J k,k' — 



«k 1 - i v k f 2 ) T (u k ' 1 + % «k' t 2 ) 



(8) 



where k is the momentum. 



Equation (JSJ shows that the scattering matrix elements 
of T k k , of the SC quasiparticles are determined by com- 
binations of the Bogoliubov coefficients Uk and v^. Note 
that the sign of the SC order parameter Ak enters the 
coherence factors according to Eq. so the momentum 
dependence of quasiparticle scattering is sensitive to the 
way the phase of Ak changes in momentum space. 
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t-matrix 


Scatterer 


Coherence factor 


Enhances 


Time reversal 


T 3 


Potential 




+- 


even 


&3 


Magnetic 


UkMk' + Wk' 


++ 


odd 


1 X todd(oj) 


Resonant 


MkHk' + Uk«k' 


++ 


odd 



TABLE I: Relation between coherence factors and time reversal symmetry properties for a set of particle conserving scatterers. 
In resonant scattering, the t-matrix is an odd-function of frequency t dd(u) — — t 0£ jd(— oj). 



Now the action of the T2 matrix on the t-matrix is 
closely related to its transformation under time reversal. 
The time-reversed t-matrix T k k / is determined from its 
transpose in the following way: 

fiv{u)=a 2 f\,^{-u)a 2 , (9) 

where for generality, we have included the frequency de- 
pendence of the t-matrix. Now since the Nambu spinors 
satisfy the relation = (*L k ) T = foh^k, it fol- 

lows from Eq. (2) that an arbitrary t-matrix satisfies the 
"CPT" relation 

T2 02 T -k',-k(~ u ) °2 T~2 = -T k ,k'M- 

Combining these two relations, it follows that the time- 
reversed t-matrix is given by 

so that if the t-matrix has a well-defined parity 8 under 
time-reversal T kk / = fi'Tk.k', it follows that 

0Tk,k'(w) = -T2 T k)k /(w) f 2 . 

With this relationship, we can rewrite the quasiparticle 
t-matrix in the form 

T k,k' = i u kUw - OvkVk')Tk,k' ^ 
+ {vkUk> + Ovk'U k )(-iT 2 )T k ^ . 

It is the diagonal elements of this matrix that determine 
quasiparticle scattering. If we restrict our attention to 



scattering processes that preserve particle number, then 
the second term in this expression is off-diagonal in the 
quasiparticle basis, and does not contribute to low-energy 
quasiparticle scattering. In this case, the quasiparticle 
scattering matrix is determined by the diagonal elements 
of 

Tk,k' = {ukUk' - OvkVk>)Tk,k' ■ (11) 
(particle-conserving scattering) 

Near the Fermi surface, the magnitudes of the Bogoliubov 
coefficients are equal \uk F \ = |«k F | = so that near 
the Fermi surface, 

Tg* oc (1 - 0sgn(A k A k ')) (12) 

In other words, time-reverse even scattering (8 = 1), 
gives rise to sign-reversing scattering, while time-reverse 
odd scattering (8 = —1) gives rise to sign-preserving scat- 
tering. 

Cases of particular interest are potential, magnetic and 
resonant scattering, respectively denoted by t-matrices 
of the form f 3 , <t 3 , and i x t odd (uj), ( where t odd (uj) = 
—t odd (~uj) denotes an odd function of frequency). The 
coherence factors and time-reversal parities correspond- 
ing to these scattering mechanisms are listed in Table HI 
The application of a magnetic field selectively enhances 
time-reversal odd scattering. We now describe a theoret- 
ical model for the field-dependent scattering of quasipar- 
ticlcs in an s± superconductor. 



B. BCS model of a Pauli-limited s± superconductor 

The Fermi surface of iron-based superconductors con- 
sists of two-dimensional hole and electron pockets cen- 
tered at r and M points, respectively^ A schematic pic- 
ture of the disconnected Fermi surface is shown in Fig. [T] 
Since the hole and the electron pockets are similar in 
both shape and volume, the interband nesting between 
these pockets may generate spin fluctuations at momen- 
tum Q = (tt, 0) which can be detected as a coherence 



peak at Q in the Fourier-transformed local density of 
states. If such a spin fluctuation induces the electron 
pairing, a sign reversal of the SC-gap function occurs be- 
tween the hole and electron pockets, resulting in s±-wave 
symmetryj^i 

We employ a simple two-band tight binding model with 
a singlet pairing. In our model, we assume an extreme 
Pauli-limited superconductor, in which the dominant ef- 
fects of the applied magnetic field enter through the Zee- 
man splitting of the conduction sea. In our simplified 
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FIG. 1: Schematic picture of the disconnected Fermi surface 
of an iron-based superconductor. For simplicity, two hole 
cylinders and two electron cylinders around F and M points 
are represented as circles. Signs of the SC gap expected for 
s±-wave symmetry are denoted as '+' and '-' in the picture. 
Arrows denote inter-Fermi-pocket scatterings. The wavevec- 
tors qi and q3 are "sign-preserving" , connecting pockets with 
the same superconducting phase, while q2 is "sign-reversing", 
connecting pockets with opposite superconducting phase. 



treatment, we assume that these effects are strong enough 
to permit us to neglect orbital effects of the field. 
The model Hamiltonian is then 



r c k,a 



£ k„ 



4,<A,, 



i Ak 1 1 



k,t -k.J, 



u k,t a -k,4. 



h.C 



(13) 



where 'c' and 'd' denote the two electron bands describing 
an electron pocket at M point and a hole pocket at T 
point, respectively. We model the Fermi surface with the 
following dispersion 

e kcr = 2t [cos(k x + k v ) + cos(k x - k y )] + p c - ap B B 
(electron pocket at M point), (14) 

£k jCT = 2i (cosfcz + cosfcy) - p d - (Tfi B B (15) 
(hole pocket at T point), 

where /is is the Bohr Magneton. Note that here, all 
orbital effects of the magnetic field on the dispersion and 
gap function have been neglected. Finally, we assume a 
superconducting field of s±-wave symmetry described by 
the gap function 

Ak = A [cos(k x + k y ) + cos(k. x — k y )] 

= 2 A cos k x cos k y . (16) 

In the numerical work reported here, we took values 
p c = 3.56£, pd — 3.2i, and A = 0.2i. With these pa- 
rameters, the minimum value quasiparticle gap on the 
Fermi surface is given by A = 1.6 A = 0.32£. 



III. TUNNELING CONDUCTANCE 

Following the ideas of Ref.— in this section we nu- 
merically model the conductance ratio measured by an 
STM experiment^. First, we review the relation between 
tunneling conductance and the local density of states 
(LDOS). Within a simplified model of the tunneling pro- 
cess, the differential tunneling conductance dI/dV(r, V) 
at a location r and Voltage V is given by 



^(r,F)oc|M(r)| 2 p(r,eV-), 



(17) 



where p(r, w) = A(r, r, ui) is the single-particle density of 
states at energy u and M(r) is the spatially dependent 
tunneling-matrix element, which includes contributions 
of the sample wave function around the tip. Thus, the 
tunneling conductance measures the thermally smeared 
LDOS of the sample at the position r of the tip. Note 
that a spatial dependence of the LDOS often arises from 
a scattering off impurities distributed in the sample. 

To filter out the spatial variations in the tunneling- 
matrix elements M(r), the conductance ratio is taken as 



Z(v,V) 



dV \ l ' 



+V) Po (eV) + 5p(r,eV) 



(r,-V) p (-eV) + 5p(r,-eV)' 



(18) 



The first part po(±eV) of the LDOS describes the 
sample-averaged tunneling density of states at bias volt- 
age ±V. The spatial fluctuations of p(r, ±eV) are given 
by the second part 6p(r, ±eV). For small fluctuations 
Sp(r, ±eV) <C po(±eV) the Fourier transform of expres- 
sion (nu is 



Z(q,V) = Z (V)(2ir) 2 5(q) 
Xq,eV) 



Z (V) 



6p(q, -eV) 



Po(eV) Po(-eV) 



(19) 



where Z (V) = 



Po(eV) 



The Fourier transformed con- 



Po(-eV)- 

ductance ratio Z(q, V) consists of a single delta function 
term at q = and a backround, which describes the in- 
terference patterns produced by quasiparticle scattering 
off impurities. The condition for small fluctuations is sat- 
isfied in the clean limit at finite and sufficiently large bias 
voltages |V| > fl- 
it proves useful to write the fluctuations in the con- 
ductance ratio as a sum of two terms, even and odd in 
the bias voltage 



Z(q,V)\^ = Z (V) \Sp+(q,eV) 



1 



1 



Sp-(q,eV) 



1 



po(eV) po(-eV) 
1 



Po(eV) po(~eV) 



(20) 

where Sp ± (q,uj) = [Sp(q, +u))±Sp(q, — u>)]/2. Depending 
on the particle-hole symmetry properties of the sample- 
averaged tunneling density of states po(V), one of these 
terms can dominate. 
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We calculate the fluctuations of the LDOS in Eq. ([20)1 
using the Green's function of a superconductor in the 
presence of impurities. As is discussed in Sec. II, the 
electron field inside a superconductor can be described 
by a four-component vector which is written as 

vpt(r,r) = ^{(r,r),^(r,r),^(r,T),-V«t(r,r)) (21) 

in real space, where r denotes the position vector and r is 
imaginary time. The matrix Green's function is defined 
as the time-ordered average, 

G Q/3 (r',r;r) = - <T T * Q (r', r)*J(r, 0)). (22) 

For an electronic system in the presence of impurities, 
the Green's function (|22[) is usually calculated using the 
t-matrix methodi 24 i 25 As is already discussed above, for a 
single impurity with scattering potential U , the t-matrix 
exactly accounts for multiple scattering off that impurity. 
In momentum space, the t-matrix is determined by the 
following self-consistent equation: 

T(k,k',w) = f>(k, k') 

+ ^?7(k,k")Go(k", W )T(k",k', w ), ( 23 ) 



where t/(k, k') is the scattering potential of a single im- 
purity and Go(k, ui) is the bare Green's function of the 
BCS superconductor without impurities. The Green's 
function for an electron with a normal-state dispersion 
£k and gap function A k in an external magnetic field B 
is 



G (k,w) = [wi-£ k f 3 -Ba 3 - A k fi] 1 , 



(24) 



Using Go(k, uj) and Eq. (|2"5|) . the Fourier transformed 
Green's function can be written as 

G(k,k» = G (k,io) + G (k,w)T (k,k»G (k». 

(25) 

Note that since translational invariance is broken by im- 
purities, the Green's function depends on two momenta, 
k and k'. 

The LDOS is determined by the analytic continuation 
G(r', r; iui n ) — > G(r', r; z) of the Matsubara Green's func- 
tion G(r', r; iu n ) = G(r', r; r)e l ( 2 "+ 1 ) 7rTT dr, 



p[T,LO) = -ImTr 

7T Z 



G(r,r;o; — iS) 



(26) 



Using Eqs. (|26|) and (|25|) . the Fourier transformed fluc- 
tuations <5p ± (q, uj) in Eq. (|20|) can be written in the com- 
pact formal 




Go(k_,z)f(k_,k+,z)G (k + ,z)], 

(27) 



where k± = k ± ^ and z = uj — iS. The q dependent 
fluctuation of the LDOS is a probe of all scattering pro- 
cesses combining two momenta k + and k with a fixed 
difference q. 

The only place that a magnetic field enters into the 
theory, is inside the bare Green- functions, Eq. (IM1) : the 
effect of the Zeeman splitting on the Green's functions 
and the t-matrices can be simply understood as a result 
of making the substitution 



a 3 p, B B. 



(28) 



This substitution will, in general, modify both the 
Green's functions and the t-matrix. Therefore, the break- 
ing of time reversal symmetry can affect the system in 
two ways: (i) a polarization of the Fermi surface, i.e., 
a different chemical potential for spin up and down elec- 
trons and (ii) for an energy dependent t-matrix (resonant 
scatterer), a magnetic field will lead to a difference in the 
scattering amplitudes of spin up and down electrons. 



IV. RESULTS 

We numerically evaluated 5p ± (q, uj) from Eq. (|2"T|) us- 
ing the model (fT3|) , considering various different types of 
impurities. For simplicity, we assumed that the scatter- 
ing is equal for the two types of electrons 'c' and 'd'. 

We computed the Fourier-transformed conductance ra- 
tio Z(q, U)| q ^o from Eq. ([T^]) by evaluating numerically 
the fluctuations <5p ± (q, ui) from Eq. (|27|i for potential, 
magnetic, and resonant scattering. The different types 
of scattering were modelled using t-matrices of the form 
listed in Table [I] The computed influence of a magnetic 
field is illustrated in Figs. [20 where Z(q, V) is plotted 
as a function of momentum q for two different values of 
the magnetic field B at a fixed bias voltage eV = A/2, 
where A is the magnitude of the gap on the electron and 
hole Fermi surface pockets. (A = 0.32t in our numerical 
calculations.) The two chosen values of B are hbB = A 
(upper panels) and pbB = A/2 (lower panels). The volt- 
age V has been chosen to coincide with the bottom of the 
quasiparticlc gap at the lower field value. 

In each of Figs. [2]|4l the value of Z is normalized with 
respect to its maximum. Furthermore, for a numerical 
evaluation of the sum over k in Eq. (|27|). we discretized 
the Brillouin zone in N = 10000 k points. 



A. Potential scattering 

At first, let us consider the case where the conduction 
electrons and non-magnetic impurity atoms interact each 
other via a Coulomb potential. Assuming the Coulomb 
interaction is screened at length scales comparable to the 
lattice spacing, we consider a local scattering potential. 
For a simple modeling of the scattering, we apply Born 
approximation, which is equivalent to taking only the 
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first term in Eq. (|23[) . The t-matrix used in the calcula- 
tions is 

f P =T 3 . (29) 

Figure [2] shows the conductance ratio as a function of 
wave vector q. We obtain a large quasiparticle inter- 
ference at the q-vectors q = (0,±7r) and q = (±7r,0) 
(denoted by q2 in Fig. [IJ connecting Fermi pockets with 
opposite sign of SC gap (sign- reversing scattering). The 
signals for the sign-preserving q-vectors q = (±7T, ±7t) 
and (0,0) (denoted by qi and q 3 in Fig. [T|) are sup- 
pressed. The obtained sign change of the SC gap is a 
characteristic feature of potential scattering which can 
also be derived from its coherence factor (ukWk' — fk^k') 
(compare Table H]). 

Potential scattering 




FIG. 2: Normalized conductance ratio Z(q, V) | q^o pro- 
duced by potential scattering, evaluated at fixed bias voltage 
eV = A/2, at "high field" /isB = A (upper panel) and "low 
field" /j,bB = A/2 (lower panel). Potential scattering induces 
predominantly sign-reversing scattering (q = q%), which is 
suppressed by a field. 

As is seen from the upper panel of Fig. [3J a polariza- 
tion of the Fermi surface by applying a magnetic field 
does not change the result qualitatively. The intensity of 
the sign-reversing peaks at q2 becomes somewhat smaller 
with increasing B. Note that although the Fermi surface 



is polarized by the magnetic field, a Coulomb potential 
is not sensitive with respect to the spin of scattered elec- 
trons. 



B. Magnetic scattering 

In addition to electrostatic interactions, if the impu- 
rity atom has a magnetic moment, there is an exchange 
interaction between the local spin on the impurity site 
and the spin of conduction electrons. For the scattering 
model, wc apply again Born approximation and assume 
a local exchange function. The corresponding i-matrix 
that is used in the numerics is 

f M = ft. (30) 

Figure [3] presents the conductance ratio for pure mag- 

Magnetic scattering 




-2-1 1 2 

q x /ji 

FIG. 3: Normalized conductance ratio Z(q, V)| q /o pro- 
duced by magnetic scattering, evaluated at fixed bias voltage 
eV — A/2, at "high field" (j,bB = A (upper panel) and "low 
field" (j,bB = A/2 (lower panel). Magnetic scattering induces 
predominantly sign-preserving scattering (q = qi.s), which is 
enhanced by a field. Note that in Born approximation, there 
is no contribution from magnetic scatterers in the zero-field 
case. 

netic scatterers. The quasiparticle interference is en- 
hanced for the sign-preserving scattering processes at q 
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vectors qi and q3 and the intensity of the peaks is al- 
most equal. On the other hand, sign-reversing processes 
are suppressed. Furthermore, we obtain a characteristic 
increase of the peak intensity with increasing values of 
B. These results are consistent with other theoretical 
works that are based on more realistic band model s 19 1 21 . 
In the zero-field case, B = LDOS fluctuations coming 
from magnetic scatterers are suppressed completely. The 
reason is the proportionality of the scattering potential to 
the electron spin. Therefore, contributions coming from 
spin up and down cancel each other when the Hamil- 
tonian is symmetric under time reversal. For B > 0, 
where the time reversal symmetry is broken and spin 
up and down electrons have different chemical poten- 
tial, coherence factors (itk«k' + I'ki'k' ) associated with a 
t-matrix proportional to 03 lead to sign-preserving scat- 
tering. Thus, for magnetic scatterers, the effect of a B 
field is driven by a Zeeman-splitting of the bulk super- 
conductor. 



C. Resonant scattering 

Resonant scattering is a multiple scattering process, 
where the main contributions to the i-matrix derive from 
higher order terms in which the energy dependence is 
determined by the bare Green's function Go(k, ui). The 
effect of a magnetic field is to replace 



Resonant scattering 
high field 




1 



q x /;t 




u — > UJ - a 3 fj, B B, 

inside all Green's functions, and thus a resonant ^-matrix 
in an external magnetic field B can be related to the 
corresponding zero-field form, as follows 



T r [oj; B] = T r [ujI - n B Ba 3 }. 
An expansion to linear order in the field 
f R (u; B) = T R {u) 1 - T' R {u) B <j 3 



(31) 



(32) 



immediately shows that the leading correction to the 
zero-field scattering plays the same role as a magnetic 
impurity. For a phenomenological treatment of resonant 
scattering, we have taken the form 



T R (u;B) 



1 



(u —iT)l- HbBuz 
(oj - iT)l + hbB& 3 

(w - 4 r)2 - (/i b b) 2 



(33) 



where we chose the width V to be equal to the supercon- 
ducting gap A at the Fermi-surface. 

The result for Z(q, V)| q ^o is shown in Fig. |4] In low 
field case fi R B = A/2 (lower panel), we obtain sign- 
preserving peaks at q = (0, 0) and q = (±tt, ±7t) of rather 
weak intensity. Moreover, some slight structures around 
these peaks appear in the QPI. Increasing the field leads 
to an enhancement of sign-preserving scattering which 



FIG. 4: Normalized conductance ratio Z(q, V)| q /o for reso- 
nant scattering, evaluated at fixed bias voltage eV = A/2, 
at "high field" (j,bB = A (upper panel) and "low field" 
HbB = A/2 (lower panel). Increasing the field leads to an 
enhancement of sign-preserving scattering (q = qi.,3) 



can be seen clearly in the upper panel of Fig. [4] where 
the intensity of the sign-preserving peaks at q = qi,3 
is increased noticeably as compared to the lower field 
case. Moreover, sign-reversing contributions to the QPI 
become strongly suppressed for higher fields so that we 
obtain a qualitatively similar picture as for pure mag- 
netic scatterers (compare Fig. [3]). This behaviour is also 
expected from the £-matrix model (|33|) which predicts for 
large enough field values a dominant term proportional 
to (T3 leading to sign-preserving scattering. 



D. Discussion 

The results in Figs. show a characteristical finite 
width of Aq ps 0.2-7T for all QPI peaks. This property 
is determined by the nesting for scattering processes be- 
tween the different Fermi surfaces (compare Fig.QJ. With 
Refi22 in mind we note that the sign preserving QPI 
peaks at q3 might be very sharp, since the nesting con- 
dition is perfectly fulfilled for intrapocket scattering. 

In the present work, we have neglected all orbital cf- 
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fects of the magnetic field on the quasiparticle interfer- 
ence. In order to justify this approximation, we recalcu- 
lated the conductance ratio for all kinds of scattering un- 
der the assumption that the magnetic field is now coupled 
to the degenerate d xz and d yz orbitals within a simpli- 
fied two orbital model for iron-based superconductors^ 
As a main result, we found that a magnetic field of the 
order of A which is coupled to the orbitals leads to a 
slight deformation of the Fermi surface giving rise to a 
small change of the nesting conditions. In all cases of 
scattering, the orbital coupling affected rather the shape 
than the intensities of the quasiparticle interference peaks 
which justifies the assumption that the orbital effect in 
Pauli-limited superconductors is very small in comparism 
to the field effect based on Zeeman splitting. 

V. CONCLUSIONS 

To summarize, we have modelled the effect of Zee- 
man splitting on quasiparticle interference in iron-based 
superconductors using an s±-wave-symmetry of the su- 
perconducting order parameter. Our model is applica- 
ble under the assumption that these superconductors are 
Pauli limited. We investigated three cases of scattering: 
The effect of a spin-split spectrum on scattering off (i) 



nonmagnetic potential scattcrers, (ii) magnetic impuri- 
ties, and (iii) nonmagnetic resonant scatterers. While 
the field has almost no effect on scattering off potential 
impurities, in the cases (ii) and (iii), the main effect of 
the field is to enhance time reversal odd scattering, en- 
hancing the sign-preserving points and depressing sign- 
reversing points in the QPI pattern. A possible scatter- 
ing model that would be in accordance with STM ex- 
periments recently performed on the chalcogcnidc super- 
conductor Fe(Se,Te) consists of potential impurities and 
resonant scatterers. However, so far there is no direct 
experimental evidence for the presence of resonant scat- 
terers in these fully gapped superconductors. Possible 
other effects that arise from breaking the time reversal 
symmetry such as the influence of a supercurrent on non- 
magnetic scatterers is an interesting open question that 
deserves future investigation. 
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